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ABSTRACT 

The diffusion of cosmic ray particles in a finite volume of space 
with a simultaneous diffusion and/or transport in energy is considered. 
The solution of the appropriate differential equation may be expressed 
as an expansion in eigenfunctions of the differential operator. If one 
approximates the solution by keeping only the lowest or fundamental 
eigenfunction one obtains the common "leakage lifetime" approximation. 
In some situations this approximation can be justified out in others 
(e.g, synchrotron or inverse Compton losses) it cannot. The reason for 
the failure in this case can be seen from the point of view of the 
expansion. The solution of the general case of Fermi acceleration, 
synchrotron losses, and energy fluctuation acting together is also ob- 
tained by this method. 
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INTRODUCTION 


There has been recent discussion in the literature as to the correct 
method of treating the loss of particles from a region of space where 
spatial diffusion and energy transport and/or diffusion are occurring 
simultaneously a common method of treating this situation when it has 
s.risen in the field of cosmic physics has been to describe it by an in- 
homogeneous, partial differential equation 


j p(£ t j ) 

J t- 



(i) 


In equation 1 is a differential operator in energy that describes 

the various energy changing processes at work within the region, is 

the average lifetime of a particle against a variety of loss mechanisms 

including leakage from the boundary, and q(E, t) describes the energy 

distribution of the particles when they are introduced into the region; 

the inhomogeneous term q is often referred to as the injection spectrum. 

Solutions of this equation are usually sought for the steady state 

case, t •= O 9 for a variety of energy transport mechanisms and 

1 2 

injection spectra. In his now classic papers * Fermi in essence solved 
this equation for the case ^ ) • In his 

first paper 1 he considered "T^r^the lifetime against nuclear collisions 

2 

of the cosmic ray particles. In his second paper , he had come to the 
opinion that diffusive leakage from the galaxy was the most significant 
loss mechanism and hence considered r=?s or the "leakage lifetime". At 
present it is not believed to be very likely that Fermi's mechanism offers 
the correct explanation of cosmic rays, however, it is generally believed 
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that he presented a correct treatment of a plausible process that should, 
in fact, occur even though it might not produce cosmic rays. 

The time independent form of equation 1 has been used extensively 
to calculate equilibrium spectra for a wide variety of problems in cosmic 
physics. Energy loss as well as energy gain has been included in and 
many different injection spectra have been considered. The use of this 
approach has been too widespread to give any references that would be 
complete, however, recently a particular application has been made to the 


case of cosmic-ray electrons 


3-12 


where it is assumed tha + the injection 

_ _/> 

spectrum is of the form - K £ and that synchrotron radiation 


and inverse Compton scattering are the main contributers to the energy 
transport term, i.e., ~ ^ ) 

predicts a solution of the form 


. This analysis 
^ for £ << 


and D(g) c < /r for g ?? ^ C. 


-/ 


where = { h> 7^i) 

The significance of this analysis lies in the fact that an observation of 

the "break" in the spectrum at E » E c (although the term "break" should be 

used with great caution as I shall show in the next section) will yield 

the product h with the result that an assumption about the value 

of b (the strength of the magnetic field or the energy density of radiation) 

will give the properties of the diffusing region through the relation 

, where R is the characteristic dimension of the diffusing 

region and D is the diffusion coefficient. 

In all of the treatments of cosmic ray electrons the inhomogeneous term 

13 i4 

is included but in some treatments of cosmic radio sources * and 
15 

X-ray sources the homogeneous equation is used. It is not appropriate 

for cosmic ray electrons for as Kardashev^ points out one can not obtain 

steady flat spectra in the presence of inverse Compton or synchrotron losses. 

14 

Melrose asserts the contrary and obtains solutions of the homogeneous 
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equation that resemble the cosmic ray electron spectrum. There appear 
to be serious errors 1 ^ in Melrose*s paper, however, which invalidate 
his conclusion so I shall consider only equations with the injection 
spectrum included hereafter. 

In spite of its wide spread use equation 1 has been recently questioned 
17 

by Jokipii and Meyer ' who quite correctly point out Jiat the analogy 
between a collision lifetime and a leakage lifetime is not a valid one. 

The probability of loss by a catastrophic event such as a nuclear 
collision may properly be characterized by a uniform probability per unit 
time but the probability that a particle will be lost by leakage at the 
boundary depends on the position of observation and even then can not be 
characterized by a uniform probability per unit time. 

At any point in the diffusing region the "ages" (time since injection) 
of the particles will be distributed exponentially with a mean age 
if collision is the dominant loss mechanism. If diffusive loss is dominant, 
on the other hand, the age distribution can not, in general, be considered 
an exponential even with a variable mean age 7?— • 

Jokipii and Meyer point out that the correct equation to solve is not (l) 
but rather 


£ f>(E, r J - y-(D ?/><&, o) /- ffa c)/l£ =%(&,£) (2) 
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where D ic the diffusion coefficient and once again we consider the 
equilibrium case, = o . The solution of this equation 

with the appropriate boundary ccnditions will yield the correct energy 
spectrum at any point. 

Earlier Shen 11 had pointed out that the diffusion term was important 

at high energies since here the leakage lifetime could be much longer than 

the lifetime against synchrotron and inverse Compton losses and the 

spectrum could therefore depend on the spatial distribution of the sources. 

Nevertheless Shen continued to treat the boundary conditions by means of 

the leakage lifetime approximation. 

Jokipii and Meyer consider a flat disk source of electrons embedded 

in a diffusing medium of infinite extent. They find that rather than 

one "break” of one power in the exponent they obtain two "breaks" of 

one -half power each at energies Zr", •- D/b ft, Ax. 

where R x and R 2 are the diameter thickness and thickness of the disk 

respectively. A similar calculation had been performed earlier by 
18 

Syrovat6kii but his results were not presented in such a way as to make 

comparison with the leakage lifetime approximation very easy. 

19 

More recently Dogel and Syrovatskii have considered a similarly 

shaped source distribution but instead of an infinite diffusing medium 

they consider a spherical region with free escape ac the boundary of 

radius R^ where R o ^R 1# They obtain results essentially identical to 

those of Jokipii and Meyer. Further departures from the leakage life- 

20 

time approximation are indicated in a paper by Longair and Sunyaev . 

21 

Furthermore Berkey and Shen claim and amply demonstrate thao in 
general the lifetime of a cosmic ray electron in a general diffusing 
region has little or nothing to do with the equilibrium spectrum when 
sychrotron and inverse Compton losses are important. Rather the 
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structure of the sources is the determining factor. 

Clearly the leakage lifetime approach does not work at all well for 
cosmic-ray electrons. This brings to mind the question of whether the 
theory of Fermi acceleration might also be in error. It turns out that 
this, and other, questions may be invetigated from a general point of 
view. If equation 2 can be solved by means of an eigenfunction expansion 
of the Green's function the leakage lifetime approximation may be seen in 
its natural mathematical setting and the question of its validity may be 
more easily investigated. 

In Section II the method of eigenfunction expansion will be described 
and the leakage lifetime approximation will be placed in its setting. In 
Section II Fermi's theory and the theory of cosmic ray electrons undergoing 
synchrotron and inverse Compton losses will be compared. We will see 
that Fermi 's theory turns out to be a very good approximation and also 
why the method fails in the latter case. In Section IV we shall see how 
the method of eigenfunction expansion may be applied to more general case 
where 4 is a Fokker-Planck operator of the form 

Z & P ' £?(*'*' ***')?' *-h-'^ i£r )/0 


In the above expression the term proportional to a x describes both Fermi 
acceleration and bremstrahlung losses and the term proportional to a 2 
describes synchrotron and inverse Compton losses. The second order 
term proportional to a 3 describes a statistical spreading or diffusion in 

energy space. This term was first considered in a cosmic ray setting by 

22 25 

Terletski and Luganov and later in a more general treatment by Davis . 
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Further discussion of the importance of this term may be found in work 

2k 13 l4 15 25 

by Morrison and it has been employe ', in many later works 9 , 

Although all of the processes discussed above have been treated by 
many authors to the knowledge of the present author the combination of the 
above mentioned Fokker-Planck energy operator with spatial diffusion, 
particle losses, and injection has never been treated before. It is 
therefore believed that the results of Section IV are essentially nw, 
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II. MATHEMATICAL METHOD 


The equation we wish to solve is 


If we define the total operator by 

t * 4- - V- 1> V f ^ 

we may write (2) as 


£ ) - 



(2 ) 


The solution may be written in terms of the Green's function 

pfc ~ j d r 'Jtf Q-f e, r j V- , r ') £ (e*r ') 


where 



( 4 ) 


and the Green's function is constructed to fit the proper boundary conditions. 
If the operator £ and its adjoint have "eigen-solutions” 
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( 5 ) 


fa* = ^A 

( 

//>: = * - e? 

where )f\ represents all of the parameters required to specify the 
solutions, and if the delta function may be expanded in these functions 
(i.e., they form a complete basis set) 


&£-£') £(r-]r') = \ pjs, r) ^(e‘ t r') 


(6) 


h 

where C" represents summation over discrete parameters and integration 
over continuous ones, then the Green's function may be also expanded by^°' ^ 


GGcM) = 5 ^ ,e ‘f - - ' 


( 7 ) 




In many cases the notion of "eigensolution" must be taken with a grain of 
salt since the solutions A will not fit the boundary conditions. 
Nevertheless, much of the time (and in particular in the cases we shall 
consider) an expansion of the form (6) can be defined in terms of transform 
theory and the analysis follows through. 

At this point we shall make the assumption that the problem may be 
separated into a spatial part and an energy part. That is to say that 
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the diffusion coefficient D is independent of energy and the energy 
operator JTm has coefficients that are independent of position. If 
this is the case, the solutions are separable 


t%V - -Py(E) 


and (5) becomes 


//J , V- - \ a )X e Is) - f (e 


/ F ' <0 


iy/s) , , , 

- A S,*' "1 (t)-fy(E) 


and dividing by f / gives 

-V-DV fi h ([) _ / +„(£ r) , 

— J- -A h * 


VO 


VO 


TZ 


Since ^ and £ are independent variables the first and second 
terms are functions of \T only and £ only respectively and 
hence must be equal to constants so we hare 

^7-DX7fiJc)+k*R k (c) - 


( 8 ) 


( 9 ) 
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rjy(B) ' V 4 l£) 



(io) 


and 


4 




(n) 


At this point we shall consider equation (9) only briefly. Solutions 
to (9) may be found for a variety of shapes of diffusing volumes. The 

usual boundary conditions are R*-o on the boundary; however, any 

* 

homogeneous boundary conditions may be applied. Although the differential 
operator is formally self-adjoint the system will be self-adjoint only 
for certain types of boundary conditions. So in general Rh 

Since the volume we are considering is finite, equation (9) will have 
solutions only for a discrete set of values of ^ . The smallest value 

fz can have will be of the order 


A, 


7 tp/l x ~ ircJl /l 


if the boundary is freely penetrated or A * where cK is the 

boundary transmission coefficient (««<, JI/l ) and is the 

linear dimension of the diffusing region. In either case we see 
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1 / 

that ho y% where is the time required to 

random walk a distance L. in the first case or the Lime required to 


k 


penetrate the boundary barrier in the second case. 

Since for higher modes one usually has kn ~ n o we see 

that particles remain in the fl ^ mode a time of the order /e / }1 L 

where 3T is the random walk time, not the boundary 

penetration time. 

We may now see what the essence of the "leakage time" approximation 
is. In this approach only the lowe st or fundamental mode is considered. 

To determine whether or not this is a valid approximation we must turn 
to an examination of (10), the energy equation. 

Since the domain of equation (lo) is semi-infinite, OS E < ^ , 
the concept of "boundary" conditions must be modified somewhat since 
setting a value of f(£) at the two end points is usually not relevant, 

r 

i.e., -f(o) ZZL const, is too restrictive and + C> is usually 

not restrictive enough. The physical condition that we must apply in 
place of boundary conditions is simply that the particle density be 
finite or that the function ■fte) be integrable 


oO 


I £r) I Ji £ 


( 12 ) 
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With this condition imposed on the solutions it often turns out 1 
(and it will indeed be the case in the following sections) that there are 
no eigenfunctions of equation (lo) that also fulfill condition (12). In 
this case we must abandon our notion of a solution in terms of "a complete 
set of orthonormal eigenfunctions." 
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On the other hand, one can often find solutions for the adjoint 


equation such that 

^ ^ ( ~) ^ /~ ^ c> ^ 

o 


for any function 


%<S) 


such that 


j l% ( *)lc»£z < . 

6> 


These solutions can be Tabled by a continuous parameter £ 



and 



The delta function / (£ -/£ ' 
of the integral 



c 


(£) ■= o 

may then be expressed in terms 



where the integral is over some contour in the complex jsr plane . In 
this situation one usually says that the operator 4 has a continuous 
spectrum of eigenvalues YO) and calls the function £ <£) 

"eigenfunction” even though this is not strictly true since 


an 
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To understand what it means to say that an equation like (13) is a 

spectral expansion of the delta function we must consider the usual 

case of expansions of a function space in its basis functions. 

26 

Using the Dirac notation ^ 0 we say that a function may be expanded 
in a basis function set and write 


//> = 2 /»><*io 

\r\ 

vherejti^ — U (x )the basis functions and the expansion coefficients are 
given by 





where the integration is taken over the domain of definition of the 
function space. Since the expansion is supposed to hold true for any 
function in the space we may de r . me the unit operator 2 and write 



f-- if 


; 1 - £ /"XV 


The delta function is the unit operator since 


i(x) 



J $(*-*'} ~ ) J * 


In the case of a continuous spectrum we can replace the summation by 
an integration and write 

//■>- JX 

and 

I - $(«-*') = fJ* !*><*/ 

What we now have of course is an integral transform where 



and the inverse is given by 

//> r f Js !s> J Js '$(>) U s (x) . 
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The theory of the Fourier, Laplace, Mellin and other transforms is 
of this type and such things a s the "continuum wave functions" of quantum 
mechanical scattering problems are to be understood from this point of 
view. In Section III we will be able to construct the delta function by 
comparing our results to known transforms such as the Mellin and Fourier 
transforms. In Section IV, however, we shall be forced to consider a 
non-standard transform based on confluent hypergeometric functions. This 
transform will be established in the appendix. 
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III. FERMI ACCELERATION VERSUS SYNCHROTRON AND INVERSE COMPTON LOSSES 

1 2 

Let us first consider the problem of Fermi acceleration. 9 In 
this case we have the equation 


U - ( —/) / 

it? ^ T~ 


The solution to (16) is 




! y /* -') 




vr , 


Ce »> r . 


-*/4 


It is easy to see that there is no value of V* for which “fs is 
integrable from 0 to c*3 . However, it is also easy to see that the 
relevant transform theory is that of the Mellin transform; 


^ - y/cj 




(16) 


(17) 


(18) 


(-9) 


7 pe'iJe' { 

V L ? 

for^ an integrable q(E'). We may therefore immediately 

t L~0 


write 



( 20 ) 


where ^<6 
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Remembering equition (ll), we have from equation (7) 



where 



W 


1 

cr~;*o 




For £< £ we may close the contour or the right hand 
infinite semi-circle without adding anything to the integral. There 
are no singularities in the right hand plane so 


= £> ■ 


For ^ we similarly close the contour to the left and 

pick up the pole at V -T V- tyf - 0 J obtain 
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If all of the particles are injected at a single, low energy E 


i * e • > / C ) — ^ L J & ~^-*V ) 

we will liave 




c, ^<c) e lltli7i) 

A # 


where 



Since, as we have seen in Section II, the mode lifetimes 
are monotori cly decreasing functions of (\ and hence 
the spectrum index — (f + ) is too. 


For this reason th» j higher the mode, the steeper the spectrum, and for 
some ew'i'W /£ & one can neglect all but the fundamental mode. This 
is illustrated in Fig. 1 where log P(s) versus log E is sketched for 
the first three modes. In the case of Fermi's origin theory, the next 

-r 

mode to be considered (the first harmonic) would be of the form E where 
so one may say that the leakage lifetime approach was quite 
Justified in this case. 

On the other hand, if we consider the case of electrons injected 
into the region with a power low spectrum, ^ (P, C)* and 

gc 28 

subsequently undergoing energy loss by synchrotron radiation ' 9 or 

29 

inverse Compton scattering , we have 


XJ= &->>&) - ^ 


(25) 


whose solution is 


JJ.UOJ.UI1 IO . 

(£) — Co* ±£ £ VAeJ 


The adjoint solution is 


fife) - c**st expf-Y/kfJ 


If we consider the variable y = E \ we see that the transformation 
in this case is equivalent to the Laplace transform and we have 


tie-*') = fa -Jt'etpp/kE - Y/u-] J({) ■ 
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Using the same analysis as before, we obtain for 


fh (*, * ') 


-/ 

= (h£ ') e *r[( fyz) - ( TZ LE) ‘ ] 


Continuing, we obtain for 




H/C) Pj£) 


and 




'e*pHxbsj'j j t ' 


This expression is well known; it has been pointed out many times 
that f^(£) has two distinct asymptotic forms. For e«^ c - (xy) 




This is the spf+rnm (with TZ.- 7 Z ) that appears in the 

leakage lifetime approximation and it has led many authors^’ ^ ^ 
to refer to a "break" in the spectrum at the critical energy, E c . 

This is, in fact, rather unrealistic, in Fig. 2 we see plotted the 
value uf the effective spectral index increment 


s= 

where 

/ 



Vi. 


GfC> 


X 


- 1/xJ ( X' / exf'U / *'] J*' 


and P — X • O . We can see that the spectral index does not sharply 
break, but changes from 2.5 to 3.5 very smoothly over about two decades 
of the parameter X = E/E . For comparison one should check the 
spectrum curves in figures of references 3,5, 8, 9> an( ^ 10. At this point it 
should also be pointed out that the spectrum of photons produced by these 
particles in the inverse Compton or synchrotron process will change its slope 
over about four decades, a fact that should be kept in mind when interpreting 
"breaks" in X-ray or gamma ray spectra. 

In this situation, the effect of including higher modes is somewhat 
more complicated. At first glance one might think that there would 
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always be a steepening at the critical energy of the fundamental mode 

r’ = Ar' i ) ^ " where R is the linear dimension of the diffusing 

fc 1'ob) ~ 

volume and for the moment we have neglected , assuming it to be very 

long compared to diffusion times. As a matter of fact. Just such an 

6 

argument was made in attempting to explain the shape of the galactic 

cosmic ray electron spectrum. Such a conclusion would seem to be borne 

19 

out by the calculations of Dogel and Syrovatskii who consider a 

spherical diffusing region of radius R q and an ellipsoidal source region 

of axes R x and R 2 where R 1 ^R q and Ri/R 2 = 10 2 . They obtain a spectrum 

that steepens by l/2 in the exponent at an energy of s and 

again by a factor l/2 at an energy on'/ h Rl . This would appear to 

indicate that R is a dimension that determines the first "break" in the 
o 

spectrum. This is misleading; it is, in fact, the dimension Ri, which 

in this problem is equal to R 0 , that determines the position of the 

first break. The size of the diffusing region has very little effect 

on the shape of the equilibrium spectrum, this is almost entirely 

determined by the size and shape of the source region. 

21 

This point he.s been made by Berkey and Shen and has been illustrated 
by them in a model calculation where a spherical source of radius 13 
in a concentric diffusing volume of radius R is considered to inject 

, /, -a 

electrons with a spectrum £)<>(£' . For this injection spectrum 
f. (E) simplifies and may be expressed in closed form, i.e. 


i( £ ) = ) 
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For completeness we shall consider essentially the same? case as 
Berkey and Shen but instead of limiting our consideration to an injection 
spectrum proportional to E ‘ we shall consider a general power law and 


make use of a convenient approximate form of f (E) . We shall also see 

m 

that to an observer inside the source region, the outer diffusing volume 
serves only as a "quantization volume" for the eigenfunctions and for 
this reason has relatively little effect on the solution* 

Consider a spherical diffusing region of radius R and a concentric 
source region of radius B with B ( • Due to the spherical symmetry 

of the problem, the eigenfunctions are ^ ) where 

b„, - 'tTfr/fi, fo(h -ri- (k*r) fnn(k m t) • If 

the source strength is Y) particles per unit volume, the resulting 
density will be 

P(f,r)*-. jr, Cn/Jb-r) {.(e) 1 

m =7 

where 


V_8/ *'• *(*-* )' 



and f (E) is given by equation (27). 

A good approximation to (27) is given by 



(29) 
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( 30 ) 



It is easy to see that t**i has the same asymptotic properties as f K 
and we can see in Fig. 2 that the spectral index is a similar function 

^ /V 

of E. The major difference between and 7*%, is that A. is over- 

whelmingly simpler. 

Now using instead of in (28), we may write 




-p 


>»-/ 




('a; ►AiJ 

where /j \ £ — [(f> m l)bE/Q^ • Since - j _ 

we may write IJO) as 


— (31) 


/r//? = ^ A 


»» 





n»3 


~ Cesfk'Bj 


r ) 

/k.'/ -0* he ) 


(31’) 


The onty role that R plays in the above expression is to determine how 
close the spacing is of the levels In this sense it acts only as 

a ’’quantization volume" as in quantum mechanics and does not have a 
profound effect on the system (of size^>-B) provided /? (( . In fact, if 
p«i\ there will be many levels contributing to a small variation of 
k m B or and we may approximate the sum in (31*) by an integral 



and R no longer appears in the expression at all. This holds true for 



/ 1 ^ 

•w (£) as well as for Pm (£) 


fty I ^ / txu wen a.£a iui 7 ^ i 7 . In fact, it would hold for 

any spectrum that is a function of k such that the integral or sum is not 
dominated by the lowest values of k. If it is so dominated, the 
approximation fa m ( * n ^ O is not valid and the 

result will be a function of and hence R since = hi = 7T//Z, 

The Fermi acceleration spectrum {•(£) can be seen to 


for kfl < < J. 


be of the kind that is dominated by the smallest values of k. 

If we note that C n 2 S B) 
and oscillates for kB>? i 

simply cutting it off at K - >/B and writing (for) X/n(far) 7Z, JL 
(assuming h<3 >• We obtain 


we may approximate the integral by 

- ! 


p<e>r) - ye 


JF 


J7T D 



A _ ro 

ye s ' 


3 7T O 






There are two asymptotic forms of (33) depending on the magnitude of k B. 


f°(e) ~ f /? A . p << ± 

) 7TO ’ 


w 
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-tptD 

y£. k e B >>i 

The critical energy is therefore seen to be determined by 

k s B rt 2 * y £c. — &/(?->) h * 

Expression (31) may in fact be summed exactly giving 

V/V/V 

f* (t > +[) \ e 8 S,\ h(Kcty - Cah (hg 3 ) 


r ke&) ~~ k&/2Cosh(fc,:&) ~1 fi'nhfke f | 

~ 7 ^J ) ~j J , r <8 




vfn; fl k • B CM ^' ' •*■****$ T jtf 


C &{h(h( t ) 


5 >» h (he F ^ 

h £ irTanUk e R) J ) 



r?£ . 


(54*) 


( 35 ) 


( 35 *) 
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If 


hgB« i 


we may write (35) 8^ (35*) 




5 % i\ (h>=c } 7 

77* n ( k/fft )' ) 



For ? 7 Z ve 



S] » h (hi t) 
k £ h 


exp(- A e#)J 


r<8 



.(pri) 

( k e£ -') r, , i 

,JJ— ertWr-B)]' r >e. 


We can see at once that for r <8 the leading terms are just those, 
apart from a numerical factor, obtained by the simple calculation of 
expressions C 34) and ( 3^'). The correction terms are small and do not 
affect the spectrum appreciably. If ?£ /J the third term in ( 36 ) is of 
the order unity but in this case the two critical energies, O/fp -/) i /? ** 
and are approximately equal and one still sees only one 

"break" at a critical energy that is determined by the size of the source 


region rather than the whole diffusing volume. 

In Fig. 3 the spectral index increment, has been 

plotted as a function of X = E/E^ where E c . - o/ff ->)!» 8 ■ f ls 

given by equation 35* r = l/2 B end R/B = 1, 2, 10 and 100. 

The primary effect of increasing R from B to B x 10 2 is to broaden 

the transition region to some degree. However, the transition remains in 
the vicinity of E c (B) and no partial break appears at a different 
characteristic energy E c (R) 

19 

From this we can see why the results of Dogel and Syrovatskii ' 

17 

agree so well with those of Jopipii and Meyer even though the latter 
considered the case 
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IV. SOLUTION FOR FULL FOKKSR- PLANCK ENERGY OPERATOR 


We now wont to consider the equation 




This equation may be written in the form 


( 58 ; 


£"+ 

" / _ 
r(k>t 

+ <$£')- r (C +?$>£) / - c> 

whe re 

Lz // 

J t / _< ) , a..J 


c - ; 

+ (y - ~ ^ 1 / J j 

If we 

make the 

subst itut ion 


t(e) ■ 

= £* Ufa, £) 


we obtain an equation for U ( c\ . iz 


t U + £[+*+b +£e] /j\[p<( 0 <~i) + hit + c. +(*+?) 6s] H-O ( 59 ; 


If we choose (\ such that 


p((p<il) + b> oS +- C ~o 


i .e . 







we may make the transformation 

2 * - - j A - + * ) & — & ^ + h 


- 30 - 



and obtain the equation 


u" +-( B~&) u'— A U- O 


(40; 


This is Kumnvr's equation 'whose solutions *.re the confluent hype rg<*on>> trie 
30 

functions ; 


UciAfce;*)- X, -£&) . 


(41) 


These functions are defined by the series 


trO 

, ff ( a; b; t ) = 5 _, 


n r O 


(A)„ 2_ 

(b)„ n; 


h 


(k2) 


which is convergant for all finite £ if & J } * J / C • 


where 


( A ) n = A(a+i}(a+Z) 


• ♦ • • 


• (Arh-l) 


and 


(A)o = 1 

An independent solution is 

u x = yX’,/f(n-A-e; z-b; 1) 

It may be seen by straightforward substitution that although there are 
two independent solutions to equation 40 and two values o TO(, namely 
<?(+ and the solutions of equation 38 have the property 

i Uj (*-, i ) - ~2 i ) 


- 3L - 


(43) 


so there are still only two independent solutions of equation 58. 
The adjoint equation is, 




( 44 ) 


s r = - *>A S " ' ^ 


where b — & * /$ i 

c r = y/tj 

By the same procedure as before, we obtain solutions 


■p(&) = /E / / r / (-/'Xt] &£ ) 


(45 ) 


The solutions depend on the eigenvalue through (?(+ ( )/) . The 

functional relationship of <?( on vAs rather complicated so instead of y‘ 
we will define a new variaole, s, with which we will label our solutions. 
Define: £ — J + S 


r 


/■ 


j = i(Z,/ai) - l 


$=- /(*&:•)*- ft. 


We now have 


I( £ )= ^ ft S] 


( 46 ) 


( 47 ) 


and 


4>> 


t) i±>s;ie ) 


( 48 ) 
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We must now determine whether or not the solutions we have chosen 


form a complete set (and if so, in what sense) and if not we must find 

the proper combination of solutions that will have this property. 

Making use of the asymptotic forms of the confluent hypergeometric 
30 

function 


?) ->i , * 


and 

? A <2 XP P ) I , R* * ? O 

we see at once that the adjoint solution has the property 


/ (i=) — y 


c?s ■ 


This means that 



(E ) < -fo (E) & £ 


will be defined only for a very restricted class of injection spectra. 

It turns out that the functions we want are the linear combinations 


f (£ J - ($E)^ e.Xf>(-£E) U(S) SXr ) 


U9) 


P(e) = (ii) * U(-s, Se) 


(50 > 


whe re 
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•JS 

-hZr 






(51) 


Using the Kummer transformation 


30 


e.xp(*) ( A j ft] (B -A- Bj -2 ) 

it is straightforward to verify that h^ i and 50 are indeed solutions of 

equations 38 and kk respectively. One may also verify from equation 51 

that these solutions are symmetric in s, i.e., 2 ? U(Sj 2 ) = 3 LU(-s t *). 

30 

The asymptotic forms of U(s, z) are"^ 


U(Sj*) 



2r 


i+fi-S 


i ir! —? bo 



so we see that 


tie) 


C -it. /£■/ 


•I •fits 

“* £ ,/£)-? O 

since - / ~*k( f "" ^i/Q 3) > if *? / ( <3 3 

• na lfis.il <i()- 3, /a,) 

S 4\ t £)<} f f -“ ) d ‘ j? e ** j 


then we have 


if * lfy( £ ) j ^ jjp^D and we have a well defined 

O 

transform. It is not clear that there is any physical reason to require 
*,< 3j , however, if we require this to be formally true we will see 
that our final Greenfe function may be quite trivially continued from 

^ to 3 so n0 re& l restriction is imposed. 

It is shown in the appendix that the proper form of the unit operator 
appropriate to this transform is. 


/ & ( x ) 


v ft -A 

X e 


HIT l (x ') ' J ) r(~; < ) 


/ <36 


n 


*L)u(s/)itx Sj xjJ 


- 1 00 
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where 


X - £, f - - 



To obtain the Green* 3 function, we note that 
and since y- fa + Y ™ ^ 

^ = -a 3 (s-s,)( ***/ 


y= Zifa**) 7 -/] 


where 


therefore 


f(x/) = 

w *>(* ' 




l ** 


X \ rh->")r(-i-fi-s ) /uus,t) K^-sX) a, 

' pm r&*s) ii-s,)(fts,) 

L 


(53) 


To evaluate the integral in equation 53 w e use the definition of *,he 
function U(s, x), equation 51 to expand the integrand as. 


( *» 



/ tii 


^7T( & 5 


I (X tJX rJU rJX ) As 
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where 




i Fh *) ^'V 

/ — — * 

a - <<)(* ) 


» 


j 



/ 7-y;/7f, * j 

(*•>,) (S +'/) 



rhi)i'(->- 


(x,j ?(’-,< iff*, ' 1 


£54) 



(s ) (5 1~ 5 1 ) 


where we have used the notation 




/ 


$ \ 1 + 2 $ 


> 



Me may now determine the behavior of each of the four terms as * 06 

from the fact 50 that as $ -,OC , —? const. We see 

immediately that the first two terms are dominated by the value of the ratio 
(///' } and that we may close the contour of integration in I to the 
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right or left as (X//' ) is less than or greater than one respectively. 
The opposite is true for II. The third and fourth terms, however, are 
dominated as £ — ^ by the gamma functions since 



so we see that III may always be closed to the right and IV may always be 
closed to the left regardless of the value of (*/x') . 

If we designate by and the contours that run up the imaginary 
axis and are closed on the right and left infinite semi-circles respectively, 
we may write if XLx ' 








We must now investigate the positions of the various singularities 
of the integrands in the $ plane. We first note that the gamma functions 
will contribute no singularities, within their respective contours as 
long as - (/+/})?&, ) . Furthermore U(s, x) is an entire function 

of s 3 and F(e,x) has only simple poles at $ — -'kfft+l) • We see, 
therefore, that in each of the integrals of equations 55 and 56 the only 
poles that contribute to the integrals are those of the denominator 
($-$>)( 5 +$ 1 ) • Evaluation of the residues is straightforward from this 
point on and if we adopt the notion^;, ( ) is the larger (smaller; 
of the two variables x and x' we may write the result in the simple form 



- * & P(-i-BfS ,yx< 

rf*Si) ( X ?J 




157 ) 


In terms of E we have 


A 


- F G 


F(-l ) / Ft. 


F(f t) ) U(- f ^;J (58 j 


y _ IfN 

At this point we may note that since 1/( ! +/3) *+ ( A 1 t / 7c)/*'>} >AW 

- 1- /3 t $. ? O for any value of ai/a 3 so the restriction *7, ( 3} 

may be relaxed and our Green's function remains perfectly regular. 


Th»' Green's function of equation 58 quite complicated and one 


would not expect simple spectra to result in general; however, in the 
case of an inverse* power low injection <b<£) £ one simple 

result can be obtained. If we examine the case of large E we may in- 
vestigate the asymptotic form of the Green's function. If E>£' 
we have 




ry*/ - i i~i,) 



(59) 




i */ {’(? s,) £ , 

The exponential quickly damps any contribution from values of E' much 

/ 

removed from E so little contribution is obtained from energies £? ^ 

/ 

On the other hand, if we have 


q(es‘) ~ fy-t-i-ts,) 

The resulting spectrum is approximately 


-2 


£ 

r (-!-(> +>') 

2? * rits.) < 


*0 




- \ 


y %- 


- Co»iK 


( 60 ) 


(61) 


This result is the same as for the case of synchrotron and inverse Compton 
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loesf s alone. This is indeed reasonable since if one assumes that f(E) 
is a power low the original differential equation, equation 58, shews 
that as the synchrotron or Compton loss term proportional to 

s* dominates the entire process. 

A more detailed examination of the asymptotic forms of F(s,z; and 
U(s,z) shows that the asymptotic forms are valid if we have 


S, ~(/i A ) 
* £ 


<< 1 


Inserting the expressions for / 2 > and £ we obtain the condition 

h^'/rc 

3 i, 

where X" is the total lifetime for the nth mode. This is the same condition 
as for synchrotron and inverse Compton losses alone so we may conclude 
that inclusion of energy diffusion and Fermi acceleration does not 
substantially alter the form of cosmic ray electron spectra at high 
energ ies . 
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V . SUMMAHY 


We have investigated the validity of the leakage lifetime approxi- 
mation that has had wide use in cosmic physics. We have seen that this 
approximation has its natural setting in the eigenfunction expansion of 
the differential equation describing spatial diffusion and energy trans- 
port. The spectra that are derived in the leakage lifetime approxima- 
tion are the energy spectra of the various eigenmodes of tne diffusion 
operator plus spatial boundary conditions. Each mode lias its own leakage 
lifetime with that of the nth mode being approximately r, ^ X/n J 
where K is the random walk time across the diffusing region of 
linear dimension L, Z V/TO D being the diffusion 

coefficient . 

We saw that for Fermi acceleration the lowest mode dominates and hence 
the leakage lifetime approximation is quite good. On the other hand 
for the case of inverse power law injection of electrons with subsequent 
synchrotron or inverse Compton losses all higher modes contribute signi- 
ficantly with tne fundamental lifetime X being of no particular signi- 
ficance. The resulting spectrum in this case is determined almost wholly 
by the spatial distribution of the sources. 

Finally we were able to obtain a solution to the more general 
Fokker-Planck energy operator in terras of confluent hyper geometric func- 
tions. These solutions are, in general, quite complex, however, for 
inverse power law injection they were seen to have the familiar property 
of steepening the spectrum by one power at high energy. 

All in all the eigenfunction expansion of the Green's function is 
seen to be an effective method for putting the leakage lifetime approxi- 
mation in its natural mathematical setting and thereby ma .ing it possible 
to investigate the limits of its validity. 



Appendix 


We wish to consider the transform 



ts ) r -/- a -} - . , 

r (ti ) } x 


and its "i iverse” 



- / 

<2 

47ri 





p{-/-e>-s) 


r(-3i) 


A** 

X U (Sj y) 



and discover to what extent we may say /- (X ) ^ s CX ) 

First we will rewrite the transform formula as a Lebes gue-StieliJes 

( L ~ S ) integral . 
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J Jq<s> 
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vhe re 


Q( 





r(u ) 



r 

rw) 


X 


-a*$ 


/ 



) 






h =0 


(- / - fi + s) ^ fy'j 


We note the following properties of 


** $-* r h , fttb*') 


X 

-? J 

0**/s-2*j n! 



3i Ut 7 S —? - M 




HU* 



(-!-S> - 
(toti) I 



**> y 

2 J * 
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If FU‘) is of unbounded variation and/or of unbounded range in X ' 
the proper definition of the l-s integral is 


<?(*>) 

f Ffr‘)J4 

Q(o) 



X -> o* 



where /r („') =/?(/') if /Ffr')/J 

F a (/') " •* / Ft*') / °r 

It is assumed that the limits exist. We may now define 


and 

respectively. 



pm 


Q<z) 

Oi») 


and 





it 6 - 



Since Q^^x') = Q (Sj * J we see that for pure imaginary s 

Q is real. We may therefore introduce two auxiliary functions of s, 


A 




to- 2 4 ( ‘>%to 




2 //' = Z /-/*> 9 '"’ 

r r o 


and 


where q^fs) is the 

/ * 

variation of Q(s,X ) over the set E of points A , 

and X = /, // f S) > & 

a.ni J r - t rH if Jr^O^-O . Likewise A*/ J Jr's) 2 O 

and lr -fr '/ Jr- < s ) * & (^(s) is real) 

where we have fl^vided the ordinate between -N and +N into R intervals. 

From the theory of the integral we know that 


P(-/-A i 


A+s) 


/v?o 




ri P (s) £ r SpL^S i) jy fa) 
— rc*s) 


We now note the fact that since for s pure imaginary X ' H (~Sj X ' ) 

/ '/ V- 5 /)y / y 

is bounded uniformly in y and s, /X // ^ /rf 

where M is a constant independent of X and s. Therefore we have 
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X. ' _ s 

I Q(*>,x*.) - )l - I j /' u(~ s j x ')dx' / 


-j/x UHX)\jy' <( Al( Xl - 


*•) 


Hr") I </f »£> 


where m E is the 
r 


measure of the set E^. If we now let R be large enough so that the maximum 


leparation of ordinate subdivisions fC-AL, = 6 


we can 


consider S’"* (s ) - C~~ ( $ ) 

X ' v >- 


We have 


2. (s) - = m(L' U) 3' ti] 

/V ' r r 

- ]Z \y h) \ <( r 


= *a?x- 

We see that as we let /? c* so that £ C3 the two functions 57 iW and 

<Gr<s) converge uniformly in £ to one another. We may define another 

/£ 

auxiliary function 


- kQ - 



and readily see that 


^r' s) = Z 4 '^ tS) 

r 



Sinae the functions CT7 r (s) and ^ (S) bracket both the 

functions /Y VT ($) and %j(s ) /'(„)_ We See that as ^ ^ ^ 
and £ 0 we a-/*/® ♦$ ) 






> % (f) rM 




* 5 r rf-/-A*s) 


Also from the uniform boundedness of )( 


-/-a -5 


^- 5 , * -) we note that 


i Uhl 

//' ' ~,i tS ) 




V/, f (-S) 

w J A/,r 


)/^ (x /J -Sw'*') *)] X ; 
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U(-SjX 
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U(-s,x')]J. 
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</A f \J*' 

o 

y" OcP 

= M ( f i fi*> A/ jy^'; /<v 

^ -Z 


Since we liave assumed that 



/<-* 


we have 



/V-+ o*> 






and the limit is uniform in s. We may now write 





U( 


XJ. 


i(-t- -■ +s) 0, 


!W 


rf? s) 


/j »v> 

A>«J 

£ -? O 
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The limits in the brackets are uniform in s, the remainder of the 
integrand outside of the brackets is uniformly bounded in s, hence the 
entire integrand approaches its limiting value uniformly in s so we may write 


/v/; 


ca*> 





we may write 


4Tt r 


l 0* 

< n 

j lx UdjX) &(*><** 

r J rf-ii) PfJl) * 


J 


and 


rtx) - 


h » i 

ft-* — 




Kr.A U > ■ 


Since the set may be covered by a series of intervals 


„fr= ZLlr,,;i 
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where 


/A ,• / 


is the lenght of the interval XV, 


! Tr„j - Xir.t Xj,r,r where X 2 r / and 

y 2 } rji are the upper and lower end points respectively of the L^l ) t-b 
interval. We have then 


where 


- £?( Xi, r i ) 

Since the sum over i may be at times an infinite seVies we must question 
whether its convergence (which is assumed) is also uniform in s. We have 


V ~ 5T ^ / - /27 2-v / £ ^ t,i j <£ xi /-*W 

- s ' / * > »' / ,-*»*/ 


This may be made as small as desired by manking n sufficiently large independently 
of s so we may write. 


-r 






2: 




r/ 






— 4 *' 
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We mast now evaluate integrals of the type 


/ •* 


3 


i (•/•&•$) i +s, pi* 




*' $(SjX*) Js 


( 


0+S 


. 'L S [ (x. I ^ fs '* ) f (ir) W-*/*) ^ A, 




A-js) 




~~ (t I 'Sj A s ; X') 


i 



r(~ 7i ) P(-t-0+^j 



It can be readily seen that the third (and fourth) term makes no cotribution 
to the integral for any value of X and • This is because all of 

the poles of this term are in the left (right) hand plane and the combination 
of Gamma gunctions always insures that the contour may be closed in the right 
(left) hand infinite semi-circle thereby enclosing no poles. We may therefore 
drop these terms and consider simply 
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X/^ 

We now note that as £ **o ; ^ / 1 f ~ ^ 

and £ (?< Xp,' _? <5 ^ h g ' ' Therefore the behavior on the 

infinite circle is dominated by the factors ( <*//„ ) if 

we may close the first (second) term to the left (right) and vice-versa if 
\ X r • X ■— we may close the contours etcher way but the 

contribution if the infinite semi circle is finite. 

Next we note that the poles at / t Z 5 - — make no contribution 

as long as the first and second terms are closed in opposite directions. 

To see this observe that when — //? the first term approaches 



The second term approaches the same value for / ~ ' ' ' - - , and 

since the two sets of poles would be encircled in opposite directions the 
contributions would cancel, term for term. So that we may ignore these poles 
we will always close the two terms in opposite directions even when X h 


¥ 
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. The first term 


We next turn to the poles at £ s. £ ( w • 
has the poles along the positive axis at - /3 ^nd the second 

term has the poles along the negative axis at £ - — + {3 . We see at 

once that when X ? % h the integral is zero since in this case each term 
will be closed to the side where it has no poles. 

If X^ X we close the first term to the right and the second 
to the left. Each term will give a contribution when $ Z £ f n J 
equal to 





* 1 






Adding these together we obtain for the integral 



y * F( 


This series may be evaluated simply by evaluating the case x 5 /„ 

If we close the contours in the same manner as for X 4 X * we will 

obtain the same series but with a correction to compensate for the contributions 
from the infinite circle. We obtain for the integral 



rrr 

\ A' (n* - ^ X J 


f . 4 ♦ 

— ✓ 7T ( 
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However, we may also close the contours in the opposite manner as for the 
caEe % p X n und pet no contributions from the poles. The contribution 
from the Infinite circle is of opposite sign so we obtain a value of Z & 
for the integral. The two methods of evaluation must yeild equal results 

so we have 

x * 

IfTTl ( ' '■> j - 7 77 1 ^ ^ > * 


or 



*1 

X F( J X ) 





Combining these results we see that the integral 


loO 


C p(~ 7 -/l s) / /-/ - 
J rfas )/’(*, $} 


x U(% X) QCSjX* ) d 5 




— o 


x>x 


-5 6 - 



, _ . 

— H-rrt e / 4 x h 

x 

- 7 7fL £ x = X « . 

If now instead of Q(% Xn ) we insert ( l r J i ) = r i) 

into the integral we will obtain Q if X is not in or an end point of 

X 

the interval [ , *47? I & if X is in the interval J r 4 l t an< i 

7 7T L G* if X is an end point of the interval p £ 

We no*: note that the intervals ][* ^ £ do not overlap and must cover 

the entire interval O to JC . Therefore if we assume that X is in 
the interval Q to the point X will fall in one of three categories; 

(l) it will lie in an interval JTr,c . (ii) it will be a boundary point 
of two adjacent intervals (with different valuses of r) or 

(iii) it will be a point of accumulation .>1' an infinite sequence of intervals. 

In the case (iii) the point X can not be unambigously assigned to any 
interval and hence the value of th6 integral is indeterminite . It should be 
noted, however, that the only way an infinite sequence of intervals can arise 
is for the function / - (X to have a point of infinite oscillation and at 

this point the function itself is indeterminite, eg., sin( x~Xp ) at the 
point X - 

For the case (i) we have ^ f\ ( X ) — -f h and for the case (ii) 

where \T jhd are the values of r 
as a common boundary point. 


we have £ ( fr + f r' ) 

for the two intervals that have V 
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We must now inquire as to the effect of taking the limit /? y £ c ? . 

A little reflection will show that for both cases (i) and (ii) the limit 


J i 

6 -*0 








This may or may not be equal to ^ (*) but it will be equal for 
"almost all" values of X (values taken by (x ) on sets of 

measure zero are ignored). It is also easy to see that case (i) corresponds 
to a point of continuity and case (ii) may (but not necessarily) correspond to 
a point of discontinuity of the function /*<* j (once again ignoring 
values taken on sets of measure zero). If X is a point of discontinuity 
the transform will take on a value half way between the values on either side. 

Taking the limit r->-o is now trivial since we have assumed that 
was Large enough so that X ^ j£" and increasing it further has 

no effect so 

h« %*(*) = - Jz » 

X">~> 


If 






exists so will //*• 

A' 



!>•* ///^ 


= F(*) 


/ ) ,rt 

6~+0 


si 

X-6 


and 
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We therefore have 


F(x) - Fix ) 


almost everywhere 
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Figure Captions 


Figure 1 - Log f>(E j vs. log E ter the first three modes of 

Fermi acceleration shoving how the fundamental mode { PJ — O) dominates at high 

energy. 


Figure 2 - Plot of spectral index increm. ■.-*? 

6 - - — f° V* - /»p (E/fft ) 

for the functions * — /° 

f= x'ctpC-ir-Vx] S / ^ ex r ltf ‘~' i/ *' J J '' 


and ^ for the approximate function ^ - y ' /(>+*) 
where - f J 

and p - P . 4* 


Figure 3 - Plot of the spectral index increment 


£ C - * ~ P 


vs. I^x ~ )>*,(£/£ c ) 



is given by equation 35 in the 


text, ^ =■ /-/^ s £T A /£ - 2, 10, and 100, 

£t - s & [(p-O kB*] ' f> = *.S~ 
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